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A Bose-Einstein condensation theory for any integer spin using noncommutative quantum me- 
chanics methods is considered. The effective potential is derived as a multipolar expansion in the 
non-commutativity parameter (0) and, at second order in 6, our procedure yields to the standard 
dipole-dipole interaction with 8 2 playing the role of the strength interaction parameter. The general- 
ized Gross-Pitaevskii equation containing non-local dipolar contributions is found. For 52 Cr isotopes 
= Cdd/tot becomes ~ 10 -11 cm and, thus for this value of one cannot distinguish interactions 
coming from non-commutativity or those of dynamical origin. 
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, Contrarily to fermions, a set of bosons can occupy a same quantum state as a consequence of the Bose-Einstein 
statistics. This is a purely quantum effect theorized by S. N. Bose and A. Einstein in 1924. During a very long time 
many people thought that this effect might be experimentally unrealizable until it was discovered in 1995 [l|. 

By simplicity, many research papers on Bose-Einstein condensation consider only bosons with spin-0 and an 
i— \\ effective contact interactions modeled by a Dirac's Delta function and then, long distance interactions such as the 
. dipole-dipole are not taken into account. 

Even though the dipole-dipole contributions are very small in current experiments, recently has been shown, using 
^-i, 52 Cr isotopes, that this interaction induces explicit anisotropy effects which can be experimentally measured [2, S 0, ll| ■ 
Although from microscopic dipole-dipole interactions one can construct effective potentials, still there are so many 
i— i \ open questions concerning to the derivation of this effective potential (or pseudopotential) in spite of the intense work 
in this field jf|. 

The importance of the dipole-dipole interaction can be seen by considering two bosons with spins Si and s 2 inter- 
acting at large distances through the potential 



ON 
O 



I. INTRODUCTION 



> 

in 

2>: y = Q ^ Sl .s 2 -3(s 1 .r)( S2 .f) y (1) 

(N ■ where r is the relative position vector, f = r/r and a is the interaction strength 0. The interaction (TTJ) cannot be 
eliminated even in the dilute gas limit [l| , and therefore (fTJ) should be included in the ultracold gas treatment. 

At short distances, considerations above are technically implemented by replacing the two-body interaction simply 
by a (5-function as the pseudo-potential in the case of contact interactions and, therefore, the effective theory is 
understood in the Gross-Pitaevskii approach, i. e the many-body Hamiltonian becomes 

^ " x , 



H = £ ( -^V 2 + «(r) ) + Y, W ^ - r i)' ( 2 ) 



i>j 

with W{y.i —Yj) = 7<5(r s ; — r^) where we have set m and h equals 1. By other hand, the wave function in this context 
becomes the order parameter for the bosonic system. 
Next step in this approach is to solve 

H *(ri,r 2 ,...,r N ) = E *(r x , r 2 , r N ), (3) 

where ^(ri,^, ...,rjv) is a symmetric wave function that in the Bose-Einstein phase one can write as 

N 

${T 1) r 2 ,...,r N ) = Y[ip(r i ) ) (4) 
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where -of course- we have assumed that all the bosons are in the ground state and therefore ^>(rj) are the one-particle 
ground state wave functions. The symmetric product ^ is known as the the Hartree's Ansatz. 
Using the Hartree's Ansatz one can write the energy of this many-body system as 



N f 

Egp = Y J d 3 r [r(r) (-V 2 + u 2 r 2 ) ^(r) + 7 (7V - 1) (|^(r)| 2 ) 2 ] 



(5) 



where the one-body potential u(r) has been chosen as a harmonic oscillator one as the trapping to confine the atoms 
in a ultracold atoms experiment. 

Equation (J5J) , by other hand, is a direct consequence of the variational method and the equation of motion obtained 
by minimizing ([5]) subject to the normalization condition 



in this context turn out to be 



Vr 2 



d 3 |^(r| 2 = 1, 



. 7 (7V-l)|V>(r)| 2 



ip(r) = pip(r), 



(6) 



(7) 



which is the Gross-Pitaevskii equation and the chemical potential [i was introduced as a Lagrange multiplier that 
takes into account the condition ©. 

In ([3]), as well as in all the discussion that follows it, except by the symmetry of the wave function, the spin effects 
are neglected and therefore, in this approximation ^(r) is a scalar. The purpose of the present paper will be to look at 
Bose-Einstein condensation theory from the non-commutative quantum mechanics point of view as it was formulated 
in [j| and derive the Gross-Pitaevskii equation in this context. Additionally we will present the derivation of the 
pscudopotential in this non-commutative quantum mechanics scenario. 

The paper is organized as follows; in Section II we review briefly our previous work on NCQM and magnetic 
dipole interactions, in Section III we extend the Hartee's approach to many body NCQM and the ideas are applied to 
Bose-Einstein condensation and generalized Gross-Pitaevskii equations are found. Section V contains our conclusions. 



where 9ij and By are constant 3N x 37V matrices. 



II. NONCOMMUTATIVE QUANTUM MECHANICS AND MAGNETIC DIPOLAR INTERACTIONS 

Usually noncommutative quantum mechanics of N particles in three dimensions is a model with a given Hamiltonian 
H(pi, Xj), with {i,j} = {1, 2 . . . , 3iV} and deformed canonical commutators (see e.g [9]), i.e. 

[xi,Xj] = iOi-j, (8) 
\pi,Pj] = iB l:j , (9) 
[xi,pj] = idij, (10) 

Although (|8|)- (fT0|) can be also implemented by using the Moyal product in the phase space, it is more convenient to 
use commutative variables instead noncommutative ones. Technically this last fact imply to write p and x according 
to the rules 



x l ^x i = Xi + ^-pj, (11) 
B 

Pi^Pi = P% + -^X]' ( 12 ) 

where Xi and pi are the standard variables satisfying the standard canonical commutators, i.e. 

[xuXj] = = \pi,Pj] , 

these rules, sometimes known as Bopp's shifts, are the starting point in noncommutative quantum mechanics. We 
would like to emphasize that the spin properties in the conventional Bopp's shifts does not appears and, therefore, if 
these effects are incorporated one must to modify these rules. 
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In reference [8j, for the case of one particle in three dimensions, the spin was introduced by positing the algebra 

[xi,pj] = iSij, \Pi,Pj]=0, (13) 

where {i,j} = {1,2,3} and 9 is a parameter with dimension of length. Equation (| 13[) is just a deformation of the 
Heisenberg's algebra similar to the Snyder one [l(| ■ 

Following [|| one realize that the algebra (|13[) can be explicitly realized in terms of commutative variables by means 
of the identification 

X-l ^ Xi — ~t~ OSij 

Pi ->Pi = Pi ■= -idi, (14) 

Si ^ Si — Sj, 

where Xi and pi are now canonical operators satisfying the Heisenberg's algebra and Si are spin matrices. Notice the 
matricial character of the non-commutative coordinate operators Xi. 

This simple observation implies that any noncommutative quantum mechanical system described by the dynamic 
equation 



id t Mt)) = H(p,x,sM(t)) 



\p 2 + V{x) 



W)) (15) 



can equivalently be described by the commutative Schrodinger equation 

id t ip(x, t) = H(pi, Xi + 0Si) t) , (16) 

where tp^t) is a spinor of 2s + 1-components. 

In |8l was discussed that the noncommutative version of the harmonic oscillator whose Hamiltonian is 



H = --V z 1 



l_a 1 

2 V + 2 X ' 

= -iv 2 + ix 2 +6»x.s + 6» 2 s 2 . (17) 

The factor 9s 2 can be absorbed by a renormalization of the ground state energy and, therefore, the new contribution 
due noncommutativity is the magnetic dipolar interaction 6*x . s. Of course if one include many particles effects one 
should reproduce the dipole-dipole interaction |T]), however our main purpose below will be implement effective 
interactions. 

III. APPROACH TO BOSE-EINSTEIN CONDENSATION FOR ANY SPIN 

Let us consider N atoms described by the Hamiltonian (|2|). In order to take into account the spin effects, we use 
the prescription in (|14[) in the two-body potential for the non-commutative coordinates f and f satisfying (|13p . i.e. 

W{r -r') = W s (r -f) + W e (r - f), (18) 

where W s and Wt are the contact and large distances contributions respectively. 

Let us now write explicitly both potentials up to order 9 2 . For the contact interaction we have 

W s {v-v') = 7 5(r-r' + l9(s-s')), 

= 7 / d 3 fce 4k Ar e i(,k As , 



= 7 5(r - r') + 6> 7 As • V r 5(r - r') + ^ (As • V r ) 2 5(r - r') + 0{9 3 ), (19) 



with As 
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The large distances term turn out to be 



Wt(r-v') = 



|r — r 
1 

|ArT 



(s-s')l' 
As • Ar 
lArl 2 



2Ar| 



(As 2 - 3(As ■ Af) 2 ) + 0(8 3 ) 



(20) 



where we have defined Ar = r — r' and Af = Ar/|Ar|. 

Since the spin contributions are dominant at large distances, one can neglect the purely coulombian contributions, 
thus the the potential at large distances becomes 



W e (r - r') « 6 



As • Ar 9 2 (As 2 - 3(As • Af ) 2 ) 



lArl 3 



lArl- 



(21) 



which contains the dipole and dipole-dipole contributions and so on. 

At very large distances, of course there are not contact interactions and the dominant term of the pseudo-potential 
is (|21l) which, compared with |T]) tells us 



(22) 



It is interesting to note that the identification (f2"2")) is similar to the relationship between NCQM and the Landau's 
problem, where 6 — 1/B corresponds to the magnetic length. This identification becomes exact at the lowest Landau 
level. In our problem at hand, identification (|22[) plays a similar role to the non-commutative Landau's one and here 
9 could be directly extracted from known data. Indeed, following Q one find that 

9 2 _^dd _ 48a ^ 2 
An m 

where ao is the Bohr's radius and m is the 52 Cr isotope mass. We therefore find 

9 ~ Kr n cm, 

which is the analogous of the magnetic field in the Landau problem. 



(23) 



A. Generalized Gross-Pitaevskii equations 

Following the Hartree's approach let us consider the derivation of the Gross-Pitaevskii equations; in this approach 
the functional of energy is 



E 



where H is the many body Hamiltonian © and is assumed that the N bosons are in the ground state. In this case 
E becomes 1 



E = d 



r ir)[ -iv 2 + \^ (r + 9s) 2 J ^(r) + 2 (|^(r)| 2 ) 2 



18 I ,3 



d 3 r>*(r) ip*(r') (As • V r< 5(r - r') + 6 (As • V r ) 2 5(r - r')) ip(r')ip(r) 



9 /" j3 j3 ' i*r \ i* f i\ f As • Ar 9 (As 2 - 3(As • Af ) 2 ) \ 
+ - A,l 3 rf(r)f(r') ^_-_V V LL U(r'Mr). 



(24) 



1 In what follows we have changed the normalization of the wave function as ip{ T ) ~~ * \fNip(r) 
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Using some identities, equation (j24|) can be written as 

E = (E GP + E c + E D ) + 0(9 :i ), 

The second contribution is 



E c = J d 3 r [-Vn • ^*(r)s^(r) + (Vn • s) (V^>* • s) ^(r) 



(25) 



(26) 



where n = |i/'| 2 is the density of particles. 

The energy Ed, by other hand, can be written 



E n = \\ d 3 r^(r)U DD (r)y^(r), 



(27) 



with 



U DD (v) = / dV^(r') 



As • Ar 9 ( A§2 - 3 ( As ' Ar ) 2 



lArl 3 2 



lArl 3 



V>(r') 



(28) 



which corresponds to the dipolc and dipole-dipole contributions. 

However, before to minimize (|25p . one should note that the depletion is very small and, therefore, we can leave 
out all terms containing derivatives of n, thus 

E c «0. 

Finally, using this last fact, we arrive to the following Gross-Pitaevskii equation 
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|V a + (r + + -L \yj(r)\ 2 + —Uod ( r ) 



V>(r) = /x0(r), 



(29) 



which, except for the dipole term ^f^pr, is just the Gross-Pitaevskii equation with the nonlocal interaction discussed 
in Q. 

The non-local equation (|29[) is complicated to solve. However one could try solve perturbatively this equation by 
assuming axial symmetry and to adding the non-local term as a perturbation one. In such case at 0-th order one find 
non-intercating vortices solutions, however the explicit inclusion oHJdd could induce drastic changes. 



IV. CONCLUSIONS 



In this paper we have studied the Bose-Einstein condensates for any integer spin by using non- commutative quantum 
mechanics as a calculation technique. This procedure induces effective dipole and dipole-dipole interactions which 
coincide with previous calculations [1, @. The comparison with ([1]) yields to ([22]) and, therefore, this identification is 
similar to the relationship between 9 and the magnetic length in the Landau problem which is another example where 
the non-commutative geometry is a useful calculation technique. This last fact follows from the identification between 
9 and Cdd/^ which is just the border between non-commutativity effects and dipole-dipole interactions. Only if (|22|) 
is holds, then one cannot distinguish interactions coming from non-commutativity or dynamical origin. 
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